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Abstract 

We investigate strong-coupling superfluidity in a uniform two-component gas 
of ultra-cold Fermi atoms attractively interacting via quasi-molecular bosons 
associated with a Feshbach resonance. This interaction is tunable by the 
threshold energy 2^ of the Feshbach resonance, becoming large as 2^ is de- 
creased (relative to 2eF, where ep is the Fermi energy of one component.). In 
recent work, we showed that the enhancement of this tunable pairing inter- 
action naturally leads to the BCS-BEC crossover, where the character of the 
superfluid phase transition changes from the BCS-type to a EEC of composite 
bosons consisting of preformed Cooper-pairs and Feshbach-induced molecules. 
In this paper, we extend our previous work and study both single quasi- 
particles and the collective dynamics of the superfluid phase below the phase 
transition temperature Tc, limiting ourselves to a uniform gas. We show how 
the superfluid order parameter changes from the Cooper-pair amplitude A to 
the square root of the number of condensed molecules {(pm) associated with 
the Feshbach resonance, as the threshold energy 2i/ is lowered. In the inter- 
mediate coupling regime, the superfluidity is shown to be characterized by an 
order parameter consisting of a superposition of A and (pm- We also discuss 
the Goldstone mode associated with superfluidity, and show how its character 
smoothly changes from the Anderson-Bogoliubov phonon in the BCS regime 
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to the Bogoliubov phonon in the BEC regime in the BCS-BEC crossover. 
The velocity of this Goldstone phonon mode is shown to strongly depend on 
the value of 2^. We also show that this Goldstone mode appears as a res- 
onance in the spectrum of the density-density correlation function, which is 
experimentally accessible. 



PACS numbers: 03.75.Kk, 03.75.Ss, 74.20.Mn 
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I. INTRODUCTION 



One of the most challenging topics in current physics is the realization of superfluidity in 
a trapped atomic gas composed of two Fermi hyperfine states. A considerable experimental 
effort has already been made to cool down Fermi atom gases, such as ^Li and [1-5]. 
The temperature can now be lowered to T ~ 0.27^, where the Fermi gas should be highly 
degenerate and the observation of superfluid behavior seems imminent [6]. 

As a promising mechanism of superfluidity with a high phase transition temperature Tc, 
making use of an atomic Feshbach resonance has attracted much attention [7-15]. The Fes- 
hbach resonance describes quasi-molecular bosons, which can mediate a pairing interaction 
between Fermi atoms. This pairing interaction is tunable by the threshold energy 2^ of the 
Feshbach resonance, and can become strong as 2i/ is decreased relative to twice the Fermi 
energy of the atoms. Using this strong paring interaction, one can hope to achieve a high 
value of Tc. Experimentally, the threshold energy 2i/ can be controlled by a weak applied 
magnetic field. Very recently, this tunable interaction was observed in a Fermi gas of ^'^K 
[16,17]. 

In our recent work [12-14] , we pointed out the importance of fluctuations in the Cooper- 
channel in considering a high-Tc superfluidity originating from the strong pairing interaction 
associated with the Feshbach resonance. We extended the strong-coupling theory developed 
by Nozieres and Schmitt-Rink [18-21], to include the effects of a Feshbach resonance and 
the associated quasi-molecular bosons. We showed that these particle-particle fluctuations 
strongly suppress Tc from the value expected within the simple mean-field BCS theory. In 
addition, the character of the phase transition was shown to continuously change from the 
BCS-type to a BEC of composite bosons (consisting of preformed Cooper-pairs and long- 
lived Feshbach molecules) as the threshold energy 2i/ is lowered. Our strong-coupling theory 
thus gave an upper limit of Tc = 0.518Tf for a Fermi gas in an isotropic harmonic trap 
potential and Tc = 0.218Tf for a uniform gas [12-14]. These values are simply the BEC 
transition temperatures, expressed in terms of the Fermi temperature Tp of one of the Fermi 
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components. 

In this paper, we investigate the BCS-BEC crossover in the superfluid state, extending 
our previous work [12-14] to the superfluid region below T^. Going past the previous BCS 
mean-field approximation [9,22-24], we include strong fluctuations around the BCS mean- 
field solution. We clarify how the order parameter described by the Cooper-pair amplitude 
A = U J2p{c-piCp^) in the weak-coupling BCS theory changes to the BEC order parameter 
related to the number of condensed bosons = {bq=o) in the BCS-BEC crossover. Here Cpa 
is the annihilation operator of a Fermi atom in one of two hyperfine states (cr I) and 6, 
is the annihilation operator of the boson molecule associated with the Feshbach resonance. 

In the field of trapped ultra-cold Fermi gases, a crucial issue is to determine a clear 
unambiguous signature for superfluidity [6,25-28]. Another important problem is how to 
experimentally track the system in the BCS-BEC crossover region. In this regard, the 
Goldstone collective mode is a very useful quantity because it is deeply related to the spon- 
taneous breakdown of the gauge symmetry associated with the superfluid phase transition. 
The Goldstone mode is known as the Anderson-Bogoliubov mode in the BCS state [29], while 
it is the Bogohubov phonon in the BEC phase [30]. In this paper, we discuss how these 
collective modes change from one to the other as we go through the BCS-BEC crossover. We 
show that the velocity of the Goldstone phonon V(j, strongly depends on 2i/, and thus it offers 
a way of observing the BCS-BEC crossover phenomenon by tuning the threshold energy 21/ 
in a cigar shaped trap (where the gas is fairly uniform in the axial direction) . We also show 
that the Goldstone mode appears as a resonance in the spectrum of the density- density 
correlation function. 

The present paper only considers the superfluid phase in a uniform two-component Fermi 
gas with an attractive interaction. In Ref. [14], we discussed the same model at and above 
the superfluid transition temperature for a trapped gas, using the local density approxima- 
tion (LDA). The extension of the present work to an inhomogeneous trapped gas will be 
considered in the future. As well known, the low energy collective modes in a trapped atomic 
gas (Fermi or Bose) are strongly altered by the trap potential. However, the equivalent of 
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the Goldstone phonon modes we discuss in this paper also arise in trapped two component 
Fermi gases. These low energy collective modes have been extensively discussed in the recent 
literature [31-33] for the weak-coupling BCS limit. In this paper, we discuss the physics of 
the Goldstone phonon mode as a function of the threshold energy 2i/. A similar analysis re- 
mains to be done for the collective modes of a trapped Fermi gas in the BCS-BEC crossover 
region. 

This paper is organized as follows. In Section II, we present our coupled fermion-boson 
model. We explain how to include the strong-coupling (fluctuation) effect originating from 
the Feshbach resonance in Section III. In Section IV, we consider the strong- coupling theory 
presented in Section III in terms of Gaussian fluctuations. In Section V, we consider the 
Goldstone mode. We first derive correlation functions describing Cooper-pair fiuctuations, 
as well as a renormalized boson Green's function for quasi-molecules associated with the 
Feshbach resonance, under the Hartree-Fock random phase approximation (HF-RPA). The 
Goldstone mode is then obtained from their poles. In Section VI, we discuss the BCS- 
BEC crossover below based on our numerical results. In this section, we also discuss the 
BCS-BEC crossover behavior of the order parameter and the Goldstone phonon mode, as 
a function of 2v. In Section VII, we discuss the coupling of the Goldstone mode with the 
density fluctuations in a gas of Fermi atoms. 



II. COUPLED FERMION-BOSON MODEL 

We consider a gas of Fermi atoms composed of two atomic hyperfine states, coupled 
to molecular two particle state. We describe the two hyperfine states using a pseudo-spin 
variable a (=t, [). The coupled fermion-boson model Hamiltonian is given by [7-15,22-24] 

per p,p',q 

+ EN, + 2i^]b% + 5rEKc_P+fiCp+fT + h.c.]. (2.1) 

Q P,Q 
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Here a Fermi atom and a quasi-molecular boson associated with the Feshbach resonance 
are, respectively, described by the destruction operators Cp„ and 6,. The kinetic energy of a 
Fermi atom is Sp = p^/2m, and + = /2M + 2i/ is the excitation spectrum of the b- 
molecular bosons. Here 2v represents the lowest excitation energy of 6-bosons, also referred 
to as the threshold energy of the Feshbach resonance. The last term in Eq. (2.1) describes 
the Feshbach resonance with a coupling constant gfr, which describes how a 6-molecule can 
dissociate into two Fermi atoms and how two Fermi atoms can form one 6-boson. The 
Hamiltonian in Eq. (2.1) also includes an attractive fermion-fermion interaction —U (< 0) 
arising from non-resonant processes [9]. 

Since one 6-Bose molecule consists of two Fermi atoms, the boson mass M — 2m and 
the conservation of the total number of particles N imposes the relation 



We incorporate this latter constraint into the model Hamiltonian in Eq. (2.1) using a 
chemical potential, Ti = H — iiN . The resulting grand-canonical Hamiltonian 7i has the 
same form as Eq. (2.1), except that the kinetic energies of Fermi atoms and 6-bosons are 
replaced as £p — > = £p — // and Sbq + 21/ ^ = ^Bq + 21/ — 2/x, respectively. In the latter 
replacement, the factor 2 in 2/^ reflects that one 6-boson consists of two Fermi atoms. 

In this paper, we investigate strong-coupling effects in the superfluid phase, as well as 
the Goldstone mode associated with superfluidity in the BCS-BEC crossover region. As a 
start, we consider a simple uniform Fermi gas and leave the effect of a trapping potential 
to future work. In this regard, we have shown in Refs. [12-14]. that while a trap potential 
enhances the transition temperature in the BEG regime, the qualitative behavior of Tc 
in the BCS-BEC crossover is not very different from a uniform Fermi gas. Within weak- 
coupling BCS theory, several papers have discussed collective excitations in a trapped Fermi 
gas with attractive interactions. (See, for example, Refs. [31-33].) 

When the Feshbach coupling term is absent in Eq. (2.1), the fermions and 6-bosons are 
decoupled from each other. In this limit, a BCS superfluid phase transition of Fermi atoms 





pa 



Q 
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and BEC transition of 6-bosons can occur, at different temperatures. These two superfiuid 
phases are, respectively, described by independent order parameters 

A = C/Ep(c-piCpt), ^^^^ 

0m = {bq=o)- 

On the other hand, when the Feshbach resonance term is present {g^ ^ 0), we find the 
following identity in the equilibrium state: 

= = = ([&o, n\) = {2u - 2;.)0^ + |a. (2.4) 

Eq. (2.4) gives [8,9] 

This last result shows that the BEC order parameter (f)„^ and the Cooper-pair order parame- 
ter A are no longer independent, due to the hybridization induced by the Feshbach resonance 
Qt-. Both A and are finite in the superfiuid phase, and there is a unique superfiuid phase 
transition in this coupled fermion-boson model. 

For later convenience, we define the following composite order parameter [8,9,12] 

A = A - g,(j)^. (2.6) 

We will find that A determines the excitation energy gap in the spectrum of fermion quasi- 
particles below Tc in the BCS-BEC crossover regime. 

III. STRONG-COUPLING EFFECT ON SUPERFLUIDITY 

1. Review on the strong- coupling theory for Tc 

In this section, we review the strong-coupling theory for Tc discussed in our previous 
papers [12-14]. This formulation is extended to the region below in the next subsection. 
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In previous work [12-14], we extended the strong-coupling theory developed by Nozieres 
and Schmitt-Rink [18] to the coupled fermion-boson model in Eq. (2.1). The equation for 
Tc was obtained by using the Thouless criterion, which states that the superfluid phase 
transition occurs when the particle-particle vertex function r(q, cu) describing the Cooper- 
channel develops a pole at q = a; 0. Within the f-matrix approximation in terms oi —U 
and described diagrammatically in Fig. 1(a), this equation for Tc is given by 

(3,) 

p ^Sp 

where 

is an effective pairing interaction. The last term in Eq. (3.2) describes the interaction medi- 
ated by 6-bosons associated with the Feshbach resonance. Eq. (3.1) is formally identical to 
the equation for Tc in ordinary weak-coupling BCS theory. However, the chemical potential 
II in the kinetic energy — Sp — /i oi the Fermi atoms can deviate strongly from ep as one 
approaches the BEC regime (where e-p is the bare Fermi energy of one spin component). 
This contrasts with simple BCS theory, where one finds that ~ sp. 

The chemical potential /i is determined by the equation for the total number of Fermi 
atoms, using the identity N — —dQ/d/j,. We include the effect of fluctuations in the Cooper- 
channel [the first line in Fig. 1(b)] as well as the Feshbach resonance [the second line in Fig. 
1(b)] in the thermodynamic potential Q. The resulting equation relating /i and N is 

N^Nl + 2N^--J2 e*''^"^ ln[l - [U - g^Do{q, zi/n)]n(q, . (3.3) 

Here /3 = 1/T is the inverse of temperature, A^p = '2J2pf{ip) and = Z], TiB(Csg), where 
f{e) and nB{s) represent the Fermi and Bose distribution functions, respectively. The last 
term in Eq. (3.3) describes the fluctuation contribution to N. The 6-boson Green's function 
is given by 

Do{r,Wr^) = -. (3.4) 

iT^n - ^Bq 



where the Bose Matsubara frequency is = 2mrT {n — 0, ±1,±2, • • •). n(q, is the 
correlation function of the Cooper-pair field operator F, = J2pC-p+iiCp+ip given by 

n(q,zi/0 = /^''ciTe-"-(T.{F,(r)F,t(0)}) 

Physically, n(q, describes fluctuations of Cooper-pairs in the normal phase. Go{q, iuJm) 
is a fermion thermal Green's function defined by 

Go{(\,iujm) = -. — (3.6) 

where the Fermi Matsubara frequency is Um — (2m -|- l)7rT (m = 0, ±1, ±2, • • •). 

The coupled equations (3.1) and (3.3) determine and 11 self-consistently. In calculating 
these equations, a cutoff is necessary to make the momentum summation converge. In Refs. 

2/2 

[12-14], we simply introduced a Gaussian cutoff e~^p''^=, which is also used in this paper. 

The self-consistent solution (Tc and /x) of these coupled equations is summarized in 
Fig. 2. When v ^ e^, since the chemical potential is at most /j, <, ep, the Feshbach-induced 
contribution to the pairing interaction / {2v — 2jjL) in Eq. (3.2) is small. In this regime. Fig. 
2(a) shows that the superfluid phase transition is well described by the weak-coupling BCS 
theory with a (weak) pairing interaction —U. In addition, we see that // ~ £f, as shown in 
Fig. 2(b). However, the chemical potential ji gradually deviates from e-p as the threshold 
energy 2v is lowered towards 2£y and below. In particular, one finds /i approaches v when 
1/ < 0. In this regime, the Feshbach-induced pairing interaction g^ji^v — 2//) in Eq. (3.2) 
is large, and deviates significantly from the prediction of weak-coupling BCS theory, as 
shown in Fig. 2(a). Since 2v is the lowest excitation energy of 6-bosons and their chemical 
potential is 2//, the situation 2v — 2//, realized in the limit of large negative values of v js-p^ 
is equivalent to the condition of BEG in a non- interacting Bose gas. Indeed, Fig. 2(a) shows 
that Tc corresponds precisely to the transition temperature of a free Bose gas of A^/2 atoms 
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when v/e-p < —1. These results show that the BCS-BEC crossover occurs in the region 
around zero threshold {2v = 0), at least for small values of g,-. This crossover phenomenon 
can be simply controlled by the threshold energy 2u of the Feshbach resonance. 



2. Strong- coupling theory below 



In order to formulate the analogous strong-coupling theory below Tc, we separate out the 
fluctuations of Cooper-pairs and condensed 6-bosons around their mean-field values denoted 
by A and 0m, respectively [34,35]. For this purpose, we write = (-F,=o)^5,o + and 
bq — 4'mSq,Q + Sbq^ whcrc Fg is defined before Eq. (3.5). Separating out the fluctuation 
contribution to the Hamiltonian in Eq. (2.1), we obtain 



2 



^ = |^ + Ee. + E Hi^P^s - ^n]% + y: ^Bqb% 



p 



+ ^ YIHPi + ^iPt] - T YliPiiPi-^ + P2,P2-,] ■ (3.7) 

Here we have introduced the Nambu field operator for Fermi atoms as = (c^^, c_pj and 
the corresponding 2 x 2-Pauh matrices (i = 1, 2, 3) acting on the particle-hole space [36]. 
The order parameter A is defined by Eqs. (2.3) and (2.6), which we can take to be real and 



proportional to Ti = 



without loss of generality. In Eq. (3.7) and the subsequent 
discussion, we use the generalized density operators = pi±g ± ip2±q, where 

P.« = E*l+?^i*P-?- (3-8) 



p 

We note that ps, = J2p,cr ^p+i ordinary density fluctuation operator. Similarly, 

one has 

Pi, = Ep[c_p_|iC^_|T + c^+|T^Uf i]' (39) 
P2q = ^Ep[c_p_|iS_|t - c^^atcl^+iJ, 

and hence 
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The operators pi, and p2q describe, respectively, the amphtude and phase fluctuations of 
Cooper-pair field fiuctuation operator Fq. In Eq. (3.7), the fermion-fermion interaction is 
seen to be neatly expressed as the sum of interactions between the amplitude fluctuations 
{—J Y^g PiqPi-q) and phase fluctuations {—jJ2qP2qP2-q)- The Feshbach resonance is also 
expressed as an interaction between the 6-bosons and the fluctuations described by p^. In 
Eq. (3.7), we have simply written pi,=o - (Pig=o) ^ Piq=o and 6,=o - 0m ^ &,=o- 

Within the mean-fleld approximation described by the third term in Eq. (3.7), the 
fermion thermal Green's function is conveniently discussed in terms of a 2 x 2 matrix Green's 
function 

G{p, lujm) = : = ,,2 ,p2 ■ (3-11) 



Here Ep = ^^^^ _|_ ^2 jg ^j^g energy spectrum of fermion quasi-particles below T^., which we 
shall call the BCS-Bogoliubov quasi-particle spectrum. Eq. (3.11) reduces to the mean-fleld 
BCS matrix Green's function [36] when the composite order parameter A is replaced by 
A. The off-diagonal (static) self-energy Ep = — Ati comes from the mean-fleld term Ati 
appearing in the Hamiltonian in Eq. (3.7). This self-energy corresponds to the mean-fleld 
diagrams shown in Fig. 3, 

^(a) ^-^i^T. ^-[^^(p, ium)] = -nU J2 ^ tanh ^Ep, 

E(b) = -r+£i^o(0, 0) T>[^-G(P, iu^m)] = -^+^7^ E tanh f (3.12) 

4P p,a;„ If - Ip p AEp 2 

p,u;m 2u- 2p p AEp 2 

where t± = Ti ± iT2. In Eq. (3.12), E(b) and E(c) include the pairing interaction mediated 
by a Feshbach 6-molecule described by the propagator Dq. E(a) comes from the (weak) 
non-resonant interaction —U. The matrix self-energy Ep = E(a) + ^{h) + ^(c) sums up to 
give 
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Sf = -Tit/eff Yl 7^ tanh (^E^ = -nA. (3.13) 

In the last expression, we have used the gap equation in Eq. (3.16). 
Using eq. (3.11), we see that 

A ^ UY.{c.,^c,^) = ^ E Tr[nG(p,^u;J] = ^/^ A tanh (3.14) 

The order parameters 0m and A can be obtained in terms of A by using Eqs. (2.5) and 
(2.6), to give 

A = A, 

^^"^ 1 ~ (3-15) 

0m = -7, \ rr 

The equation for the composite order parameter A = A — gr4>m can then be rewritten in the 
form 

A = t/eff E ^ tanh (3.16) 

where Ues is defined in Eq. (3.2). This self-consistent equation has the same form as the 
BCS gap equation if we replace A — > A and // — > £f- We also note that Eq. (3.16) reduces 
to the Tc equation in Eq. (3.1) when A ^ [12]. 

It is important to emphasize that the 2 x 2-matrix single-particle Green's function in 
Eq. (3.11) only includes self-energy effects arising from the (off-diagonal) static mean-fields 
produced by the Cooper-pairs and the Bose-condensed 6-molecules. In the approximation 
we use in this paper, the frequency-dependent fermion sclf-cncrgies associated with the 
order parameter collective modes are not included in Eq. (3.11). However, Eq. (3.11) does 
implicitly involve the self-consistent renormalized values of A and as determined by the 
order parameter fiuctuations (for further discussion of this kind of approximation, see Ref. 
[37]). An improved theory of the BCS-BEC crossover would be based on including the 
fermion self-energies arising from coupling to collective modes [38] . 

The chemical potential /i is determined from the equation for the total number of parti- 
cles N. As in our discussion of the strong-coupling theory for Tc [12-14], we work with the 
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thermodynamic potential ft consisting of a static mean-field part (= JImf) and a fiuctuation 
part (= SQ) originating from the particle-particle Cooper-channel, as modified by the Fesh- 
bach resonance. The self- consistent equation for N is given using the identity N — —dfl/dii. 
The mean-field part is easily obtained from the first four terms on the right hand side in 
Eq. (3.7) [39], 

^^MF = ^ + E(^P - E,) - 2TE ln[l + e-"^^] + ln[l - g-^^-"]. (3.17) 

^eff p p q 

Our approximation for the fiuctuation contribution 50, below T^., corresponding to the 
contribution at shown in Fig. 1(b), is given diagrammatically in Fig. 4(a). Among 
the interaction terms in Eq. (3.7), we first carry out perturbative expansion in terms of 
~ J J2q P2qP2-q, wfilch descrlbes the interaction between the phase fiuctuations of A. Sum- 
ming up the loop-type diagrams in Fig. 4(a), we obtain the phase fluctuation contribution 
to Sfl 

<^^2 = ^ E + ^n°2(q, lUn)] . (3.18) 

In this i/„-summation and in the following equations, we omit the important convergence fac- 
tor e"^"'^, for simphcity of notation. The generahzed density correlation function n22(q, ii^n) 
is defined by [34,35] 



^iM^lVn) = -\ dT^^-^{T,{p,^{T)p,.q{m (^,J = 1,2,3) 

1 



= ^ E TrhG^(p ^, iLOrn + ii^„)TjG(p - ^, iu;„)l , (3.19) 



2 

where the second line (= 11° ) is the approximation neglecting the effect of the interactions 
—U and g-r. 

Eq. (3.19) also defines other correlation functions, which will be important, such as Hn 
and ni2. Physically, rin and 1122 describe, respectively, the amplitude and phase fiuctuations 
of Cooper-pairs. 1133 describes density fiuctuations in the gas of Fermi atoms. Iljj with j 
describes a coupling between fluctuations, e.g., is a coupling of amplitude fluctuations 
with the phase fluctuations (amplitude-phase coupling). 
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Summing up the Matsubara frequencies in Eq. (3.19), we obtain [34,40] 



P 



Hi- 

p 
p 

p 



Vf^p- 


2 




^p+f^p-f - 




^P+f + ^p-f 


^p+f^p- 


2 


J* 1 2 "2 








^p+f^p- 


2 


{E^+i - E^_^y + 


ep+fep-f+A% 


^P+f + ^P-f 


^P+f^P- 


q ^ 

2 


(^P+f +^P-f)^ + ^^ 


^P+f ^P-f - 




^P+f -^P-f 


p+2 P- 


2 


(^P+f-^P-f)' + ^^ 




A^ 


^P+f +£^P-f 


^P+f^P- 


9 ^ 

2 





(3.20) 



1 - fiE, 



p+f. 



(3.21) 



/(^P+f)-/(^p-f). 

i-/(^p+f)-/(^p-f; 



(3.22) 



The correlation functions (i 7^ j) describing the couphng of different operators are given 

by 



n' 



- _\^f ^P+2 '^P 2 



V-/^P+f , ^P 



^p-f^(^p+?-^p-?)'+^^ 



p+2 P~2 ^ P+2 "P 2 



P ^^P+f E^_J(E^^,+E^_,)^ + ul 

~ " ?^-^P+?/2 Ep_g/2^ {Ep+q/2 - Ep-q/2)^ + Vl 



fiEp+^j-m.,) 



1 - /( V§) - HE,..] 



(3.23) 



f {Ep+q/2) - f{Ep_g/2) 



1 



^^13 



-AE 



-Ei,+g/2 -Ep-g/2 ^ {Ep+q/2 + -Sp-g/2)^ + 
^p+9/2 + Cp-g/2 -E'p+g/2 " -E'p-g/2 



[l-/(ii;,+,/2)-/(^p_,/2)], (3.24) 



-AEE 



Ep+q/2Ep-q/2 {Ep+q/2 " -Ep-g/2)^ + 
Cp+g/2 + Cp-g/2 Ep+g/2 + -E'p-g/2 



P P 



Ep+q/2Ep-q/2 {Ep+q/2 + -£^p-g/2)^ + 



f {Ep+q/2) - f{Ep_g/2) 

l-f{Ep+,/2)-f{Ep_,/2)], (3.25) 



with 



-n?9, IlqT = —1123 ^^'^ ^31 ~ ^13- The density-density correlation function 1133 



12 5 -^-^32 

and the related coupling correlation functions 1153 and will be used in Section VII, where 
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we show that the Goldstone mode describing the collective phase oscillation of Cooper-pairs 

has spectral weight in the density correlation function. This weight comes from the coupling 

to the amplitude-density (H^s) and phase-density (n.23) correlation. 

The factor E^^^ ~-^p-f denominator of the first hne in Eqs. (3.20)-(3.25) describes 

scattering between excitations with momenta p ± q/2 in the same quasi-particle band Ep. 

For this reason, the first line in Eqs. (3.20)-(3.25) is referred to as the intraband term [34]. 

Since the thermal excitations of fermion quasi-particles are absent at T = 0, the intraband 

term vanishes at T = 0. On the other hand, the second line in Eqs. (3.20)-(3.25) is finite 

even at T = 0. The factor E„,q + E„_q in the denominator describes interband scattering 

i'~<~ 2 V 2 

between -E'p+a and —E^_i, and thus the second line in Eqs. (3.20)-(3.25) is called the 
interband term. The intraband term is known to give rise to Landau damping of collective 
modes below the excitation gap 2A, while the damping due to the interband term only exists 
above 2 A [34]. We note that the gap equation in Eq. (3.16) can be neatly expressed in terms 
of in Eq. (3.21), namely 

l + ^n°2(0,0) = 0. (3.26) 

The fluctuation contribution {5VLi) involving the amplitude fluctuations of Cooper-pairs 
is similarly obtained by summing up the loop-type diagrams in Fig. 4(a) associated with 
the amplitude-amplitude interaction — j ZlgPigPi-«- This gives 

'^^i = E ln[l + ^n?,(q, iv^)] . (3.27) 

The phase and amplitude fluctuations are coupled with each other through the amplitude- 
phase coupling 1152 in Eq. (3.23). This additional coupling effect can be formally incorpo- 
rated into 5Vti by replacing n5i(q, Wn) in Eq. (3.27) by nii(q, ii/„), where 

Ii,M.i^n) = n?i(q,^z/„) + n02(q,^i/J— — (3.28) 

The second term describes the amphtude-phase couphng effect through the couphng correla- 
tion functions n'j'2 and Ilgi. Eq. (3.28) is obtained by summing up the fluctuation diagrams 
shown in Fig. 4(b). 
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In summary, the fluctuation contribution 5Qu to the thermodynamic potential involving 
only the non-resonant interaction —U is the sum of Eqs. (3.18) and (3.28) with U^^ — > flii. 
The sum is given by 



^^u = ^ Elii[[l + ^tlu{q,ti^nm + 7n°2(q,^i/n)] 

= ^ + ^n?i(q,zi/J][l + ^lil^ia^iv^)] - (^)'n?2(q,^i^n)n°i(q,zi/0 



= :4Elndet[l + ^fl°(q,^z/n 



(3.29) 



where we have used the well known identity det A — e^^^^"^^ in the last expression. n°(q, iun) 
is a 2 X 2-matrix density correlation function, defined by 



fl°(q,ii/„) = 



^ Y[0 ttO ^ 

'-hi '-'-12 



17° 17° 



(3.30) 



Finally we consider the fluctuation contributions 5J1fr from the Feshbach resonance term 
J2q[blPq +bqPg] lu Eq. (3.7). First we sum up the diagrams described in Fig. 4(a), where 



the dashed line now represents the 6-boson Green's function DQ{q, iun) in Eq. (3.4), to give 

(3.31) 



^^FR= ^ETrln l-^i?°(q,^^^n)S°(q,^^/n) 



Here i)°(q, ii/„) is a 2 x 2-matrix 6-boson Green's function defined by 

1 iz/„r3 + [,Bq 



-.2 ' 



(3.32) 



i^nTz - + Hq 

and S'^(q, iz/„) is a 2 x 2-matrix correlation function (neglecting the effect of —U and ^fr) 
defined by 



S°(q,ii/„) = - / dre 



^ Pqir)pt,iO) P,(r)p:,(0)^ 



(3.33) 



\ptir)ptq{0) P+(r)pI,(0) ; 

Using the definition = ± ip2±q, we find that S° can be expressed in terms of the 
matrix elements of n° as follows: 
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2° = 2W-^U°W 

^n?, + no, + i(n?,-nOi) 



17° - 17° 

^^11 '-'-22 



V 



'-'-U '-'-22 



(3.34) 



where Vt^ is the unitary matrix 



V2 



^ 1^ 



(3.35) 



Next we renormahze the fluctuations in Eq. (3.31) by including the effects of the non- 
resonant interaction —U on the correlation function S°(q, ivn), working within the HF-RPA. 
This results in S°(q, ii/„) in Eq. (3.31) being replaced by Eu, where [34,35] 



(3.36) 



The resulting expression for 5Qfr involves the fluctuation effects related to both the non- 
resonant interaction —U and the Feshbach resonance coupling parameter g,-, namely 



1 ^„ . g^f^o, 



<^^FR = T^ETrln l-^D°(q,zi/n)S[;(q,^z/„ 



= ^ETrln{l-fD°(q,.z.„) 



1 + jS°(q,ii/„ 



''sO(q,^z/„)}. (3.37) 



The total fluctuation contribution SQ is given by the sum of Eqs. (3.29) and (3.37). Re- 
ing the definition in Eq.(3.34) and the re 
Trln[l -I- jH°], this sum can be written as 



calling the definition in Eq.(3.34) and the relation Trln[l + |n°] = Tr ln[iy{H-^E:°}iy~^ 



= ^ E Tr ln{l + ^ [t/ - g^D^q, tuM'{cL, tr^n)}- 
Putting everything together, the total thermodynamic potential = Qmf + is 



(3.38) 



cff 



E(^P - E,) - 2r E In [l + e-^^^] + T E In [l - e 



-0a 



+ i-E^ln{l + ^[t/- ^,'^°(q,ii^n)]5°(q,i^n)}. 
The equation for the total number of particles is then obtained from N — —dfl/d/i. 



(3.39) 
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In taking the derivative with respect to /i, we note that the order parameter A also 
depends on the chemical potential /i, in addition to = Sp — fj. and ^Bq = + 2i/ — 2/1. 
Thus one needs to calculate 

dn dA _ dCtuF dA dSn dA 
dK~dpL ~ dA dn ^ 'dK'dpL 

= (3.40) 

dA di^ ^ ^ 

In the last hne, we have used the fact d^Mp/dA — 0, which holds when A satisfies the gap 
equation (3.16). However, as will be shown in Section IV, since d5Q/dA in Eq. (3.40) is a 
higher order correction within the perturbative approximation we are using, we can neglect 
the contribution from Eq. (3.40). Thus the dependence of A on // only leads to higher order 
corrections. The resulting equation for N is 



N^2<l>l + Y.[l-f tanh ^Ep] + 2 



p 

- ^E|;Trln{l + ^[C/-^.^^°(q,ii.„)]H°(q,z^„)}. (3.41) 

Here it is understood that the /^-derivative in the last term only acts on the chemical potential 
involved in ^p and ^Bq- Eqs. (3.16) and (3.41) give us the required self-consistent coupled 
equations for A and /j, in the superfiuid phase below T^. We will discuss our numerical 
self-consistent solutions of the coupled equations (3.16) and (3.41) in Section VI. 

When one neglects the fiuctuation contribution given in the last term in Eq. (3.41), 
we obtain the mean-field expression obtained in Ref. [24] in the context of high-T^ super- 
conductivity. We also mention that Eq. (3.41) reproduces Eq. (3.3) for the normal phase 
T > Tc, where A = = 0. To see this, we note that the second term on the right hand 
side in Eq. (3.41) reduces to 2A^^ in Eq. (3.3). In addition, since the phase and ampli- 
tude fluctuations are indistinguishable when A = 0, we flnd H^^^ = TI22, and thus Sq in 
Eq. (3.34) becomes diagonal at Tc. Noting that H5i(q, = Tc) = — 4H(q, ii/„) and 

E22{ci,iUn,T — Tc) — — 4H(q, —ivn), where H(q, ii/„) is given in Eq. (3.5), we find that the 
last term in Eq. (3.41) (= 5N) reproduces the last term in Eq. (3.3). More explicitly, we 
have 
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(3.42) 



where we have used Dl^{q,ih'n) = -Do(q, ^z^n) and D22(q;^'^ra) = -Do(q, —ii^n) in the last line. 
Thus Eq. (3.41) is equivalent to Eq. (3.3) at Tc. Our present strong-coupling theory giving 
the coupled equations (3.16) and (3.41) for the superfluid phase is seen to smoothly go over 
to our previous discussion at and above [12-14]. 

Each term in Eq. (3.41) has a simple physical meaning, which is useful to discuss. The 
first term 



2N^ = 2cj>l = 2(6,=o) 



(3.43) 



gives twice the number of Bose- condensed 6-bosons. The second term 



(3.44) 



describes the number of Fermi quasi-particles. This expression can be directly obtained 
from A^F = J2p,a{(^l,a^pa) Eq. (2.2). To understand the physical meanings of the last two 
terms in Eq. (3.41), it is convenient to divide the //-derivative in the last term into the 
derivative acting on ^Bq = g^/2M -|- 2i/ — 2// in the 6-boson Green's function -D°(q, it/„) and 
the derivative acting on = Sp — /j, in S°(q, ii/„) (= d/diip). Using the identity 



2 E M^Bq) = - E [l + ^ E Trp°(q, ^^n)] 



(3.45) 



we can write Eq. (3.41) as 



N^2N^ + Nf + 2N^ + 2Nc, 



(3.46) 



where and A^c are defined by 

2iV^^-^[l + i^Tr[D(q,z^„) 



2Nr 



^ E ^Tr ln{ 1 + i [t/ - g^D%q, S° (q, } 



(3.47) 
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Here D{q, ii/n) is a renormalized 2 x 2-matrix thermal 6-boson Green's function defined by 

1 



-D(q,ii^„) 



The 6-boson matrix self-energy 



(3.48) 



E(q, = j5[/(q, = j [l + ^^n] 



S°(q,ii^n), 



(3.49) 



describes the Feshbach coupling with Fermi atoms. 

The renormahzed 6-boson Green's function in Eq. (3.48) has the same form as that 
obtained by Kostyrko and Ranninger [24] calculated within the HF-RPA. Comparing iVg 
in Eq. (3.47) with Eq. (3.45), we may interpret as the number of non- condensed b- 
bosons, as renormalized by the Feshbach resonance. In analogy to our previous discussion 
of the strong-coupling theory for Tc [12-14], Nc in Eq. (3.47) may be understood as the 
fluctuation contribution to N from Cooper-pairs associated with the (dynamical) pairing 
interaction given by 



C/eff(q, iun) = U - g^D^iq, ivn) 



(3.50) 



We note that t/eflf(q — ^li^n = 0) = t/effl, where C/eSf — U -\- g^/{2v — 2/x) is the effective 
pairing interaction appearing in the gap equation in Eq. (3.16). 

The renormalized 6-boson Green's function in Eq. (3.48) can be also written in the form 

1 _ iz>„T3 + isq - «Ti 



ll^nTs - Cb, - KTi + e,Bq 

where the renormalized parameters are given by 



(3.51) 



i—\Iiui2 - 11(721], 



g, 



es, = eB,+ ^[nt;n + nt,22], (3.52) 

The correlation functions liuij = 1, 2) are the matrix elements of flu defined by 

^[/(q,^^^n) = [l + ^n°(q, ^i/n)] ~'n°(q, ^^^n)• (3.53) 
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More explicitly, we have 



n 



- YT^n^ 



22 



(3.54) 



42 



[i + fn?,][i + fnoj-^nM,' 

n[/2i(q,i^'n) -n[/i2(q,ii/„). 

Here Hn is defined in Eq. (3.28), and 1122 is similarly obtained from Eq. (3.28) by inter- 
changing 1 2. 

When we take q — iun — 0, the denominator of Eq. (3.51) reduces to 

a..„ - <o,or ^ (1 + f n»,(o,o)l[i + ^^n;,(o,o)l^^^^^|^. (3.55) 

The expression in Eq. (3.55) clearly vanishes when the gap equation in Eq. (3.26) is 
satisfied. This means that the excitation spectrum of the renormalized 6-bosons described 
by i) is always gapless at q = for T < Tc. This is a desired result, because the Bose- 
condensation of 6-bosons characterized by should be accompanied by the appearance of 
a Bogoliubov phonon (Goldstone) mode having a (uniform system) gapless dispersion. Thus 
the strong-coupling theory presented in this section correctly includes a gapless spectrum for 
the symmetry breaking Goldstone mode. This result is not obtained if we use a static mean- 
field theory, neglecting the fluctuation contribution given by the last term of Eq. (3.41). 
Within such a mean-field theory, the &-boson excitation spectrum is given by + — 2/x, 
which always has a finite excitation gap at q = 0. This gap is given by 5E = 2v — 2^ — 

2 

— ^n[/22(0, 0) > 0, assuming that the gap equation in Eq. (3.26) is satisfied. We note that 
nc/22(0, 0) is defined in Eq. (3.54), with 1122(0, 0) = lil^i^, 0) since n?2(0, 0) = 0. 

The gapless behavior of the renormalized 6-boson can be also verified formally by noting 
that the Hugenholtz- Pines theorem is verified, namely [24,41-43] 



g^r n°,(o,o) 
2 i + fno2(o,o) 



En(0, 0) - Ei2(0, 0) = u^,^L. = 2/. - 2v. (3.56) 



Here we have used the gap equation given in Eq. (3.26). 
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In order to discuss BCS-BEC crossover in the superfluid phase, one needs to solve the 
coupled equations (3.16) and (3.41) numerically. Results will be discussed in Section VI. 
Here we briefly discuss limiting cases which can be treated analytically. 

When the threshold energy 2v of the 6-boson excitation spectrum is very large (BCS 
limit), 2v ^ 2/x can be realized because the chemical potential ii is at most the order of 
e-p. In this case, the Fermi atoms are dominant particles while the effects of 6-molecules 
described by ^m, ?T'b(Cb9) ^iid Eq. (3.41) can be neglected. In addition, since the last 

term in Eq. (3.41) is small when one is dealing with a weak non-resonant interaction — [/, 
one can also drop this term. Thus we see that Eq. (3.41) reduces to = A^f (where N-p 
is defined in Eq. (3.44)), which simply gives // ~ as far as e-p » T. When this result 
is substituted into the gap equation in Eq. (3.16) with [/eg — > U and A — > A (note that 
2i/ 3> 2yLi.), one obtains the usual BCS gap equation for the Cooper-pair order parameter A. 

In the opposite hmit 2v < —2ep (BEC hmit), since the 6-boson branch has an energy lower 
than the two fermion band energy, most Fermi atoms will combine to form 6-molecules and 
hence the fermion correlation functions Iljj becomes less important. Then the gap equation 
in Eq. (3.26) can be rewritten as [see also Eq. (3.56)] 



g^r n°,(o,o) 
2 i + fno2(o,o) 



2/. ^2.+ ^ ^ (3-57) 



This says that the chemical potential has the energy of the bottom of the 6-boson excitation 
spectrum. Substituting 2// = 2v into Eq. (3.41) with S'^(q, ii/„) = 0, we obtain 

y = 0n. + E^B(|^). (3.58) 

This is just the equation for BEC in an non-interacting uniform Bose gas with N/2 bosons 
of mass M. Thus the coupled Eqs. (3.16) and (3.41) reproduce both the BCS-phase and 
BEC-phase for two limiting values of the Feshbach molecular resonance threshold 2v. 

We again remind the reader about the many body approximation our whole discussion 
is based on. As we noted after Eq. (3.16), we have included the fluctuations of the order 
parameters A and in determining self-consistently the values of N and //. However, we 
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have not included the exphcit self-energy of the fermions due to these dynamic fluctuations. 
For further discussion, see Refs. [37,38] . 



IV. RELATION TO A GAUSSIAN FLUCTUATION THEORY 



It is interesting see what type of fluctuations are taken into account in the strong-coupling 
theory presented in Section III. In this Section, we show that the coupled Eqs. (3.16) and 
(3.41) describe Gaussian fluctuations around the mean-field order parameters A and ^m- To 
see this, we employ the functional integral formalism because it treats fluctuations in a very 
clear fashion. Indeed, using this formalism, Randeria and co-workers [20,44] showed that 
the strong-coupling superconductivity theory developed by Nozieres and Schmitt-Rink [18] 
can be understood as a Gaussian approximation for the fluctuations in the Cooper-channel. 
Very recently, Milstein et al. [15] also employed this approach to the coupled fermion-boson 
model in Eq. (2.1) and reproduced the equations for Tc given in Eqs. (3.1) and (3.3). 

The partition function for the coupled fermion-boson model in Eq. (2.1) in the functional 
integral formalism is given by 



where p = V/i and x = (r,T). "ifai^) is a Grassmann variable for Fermi atoms and 
^l{x) is its conjugate. $(x) is a complex scalar Bose field describing 6-molecules. When 
we introduce the Cooper-pair field A using the Stratonovich-Hubbard transformation, the 
partition function in Eq. (4.1) is transformed to 




(4.1) 



The action S has the form 



S 





(4.3) 



23 



where the action S is now given by 

r r\A(x)\^ 



Jo J U 2 

P ^9 



(4.4) 



Here we have introduced the Nambu representation for fermions as = (^j.^j). The 
corresponding 2 x 2-matrix thermal Green's function G{x) in the (r, r)-representation is 
given by 

^2 



G-\x)^-dr-{^-^i)n + A{x), 



(4.5) 



A{x) 



^ A(x) ^ 



A(xY 



(4.6) 



where 

/ . \ / 

A{x) - gr^ix) 

A{xy - g,^^{x) 

We have also introduced a matrix notation for the 6-bosons as l>^(a;) = {^"^ {x) , ^{x)) , as 
well as the corresponding matrix 6-boson Green's function in the (r, T)-representation, 

D'{x)-' = -d^T, - + - 2/.). (4.7) 

Since the action S in Eq. (4.4) has a bi-hnear form with respect to ^ and one can 
immediately integrate out the fermion degrees of freedom using the formula 



J V^lV^^e-^^^^'^^^'^' = det[M], 



(4.8) 



where i and j represent the variables (r, r, cr) involved in ^(^(r, r). After the functional 
integration, the action is reduced to {= S) 



^= -Trln[-G-^] +^ dr j 



dr 



- -^{xyb\x)-^^{x) 



(4.9) 



The trace in Trln[— involves the integrals over r and r, as well as a trace over the 
particle-hole space. Eq. (4.9) has a very physical form because it is described by two kinds 
of Bose fields A[x) and ^{x) directly related to superfluid fluctuations. 
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In order to extract the fluctuation contribution from Eq. (4.9), we divide the Bose flelds 
A(x) and $(x) into their mean-field values (A, (f)^) and fiuctuations from these values. In 
the functional integral formalism, the mean-field approximation corresponds to the saddle 
point solution, given by 

TTT^ = , , = 0. (4.10) 

Since we are working with a uniform system (no trap), we can take a uniform (and real) 
solution as the saddle point solution. Then Eq. (4.10) gives the two equations: 

A = C/G'i2(x) = ^ E G'i2(p,iu;J = ^E^^tanh^E^p, 
(2u - 2//)0^ = -grGuix) = -| E ^ tanh ^E^, 

where Gi2{p,iu!m) is the off-diagonal component of the fermion Green's function in Eq. 
(3.11). Eq. (4.11) is equivalent to Eqs. (3.14)-(3.15). 

Writing A{x) — A + SA{x) and ^{x) — 4>ja. + S^{x), we can express the partition function 
Z in terms of the fiuctuation fields SA{x) and 5^{x) as 

Z = e-^° J VS^^VS^VSA^VSAe-^^ . (4.12) 

Here the action Sq does not involve the fiuctuation fields, being given by 

^0 = /3 E - f E U - 1^ H-Gi^r'] + 
= /5Eep-f EU- E^M-^^(P'^^-)"'] + /5^- (4-13) 

Here and in the following equations, G{x) is given by Eq. (4.5), where A{x) is replaced by 
the saddle point solution A = A — gr4>ra- In the last expression, Gijp^iujrn) is given by Eq. 
(3.11) and the trace is over the particle-hole space. (The trace in the first line involves the 
integrals over both r and r.) In Eq. (4.12), 5S describes fiuctuations around A and ^m, 



5S ^ t dr I rfrf^^^ - U^'^{x)b\x)-^5^{x) \ + V Lil^IV[(G'(x)5A(x))"], (4.14) 
Jq J ^ U 2 J n 
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where 6 A. has the same form as Eq. (4.6), in which A(x) and $(x) are now replaced by the 
corresponding fluctuation variables 6A[x) and 6^{x), respectively. We note that since we 
expand A{x) and ^{x) around the saddle point solutions (A, ^m), the first-order terms in 
5A(x) and 6^{x) vanish. 

We now employ the Gaussian fluctuation approximation to Eq. (4.14). This approxi- 
mation corresponds to only retaining quadratic terms with respect to the fluctuation flelds, 
giving 

5S^ dr [ dr\^-^^^^ -ls^Ux)b''(x)-^5^x)] + Itt\G(x)5K(x)G(x)5K(x)]. (4.15) 
Jo J ^ U 2 -I 2 L -I 

Eq. (4.15) involves d^S/dA'^ while neglects higher order terms d'^S/dAP- {n > 3). This 
point is crucial when we calculate the equation of the total number of particles from the 
identity N — —dVt/dii: The thermodynamic potential Vt = —TlnZ has the form 

Q = TSo-TlnJ V5^^V5^V5A^V5Ae-^^ 

= TSo + A{A,cl)^,l^). (4.16) 

Since the last term in Eq. (4.16) is 0{d'^S / 9A^) , the derivative of this term with respect to A 
gives a higher-order correction under the Gaussian fluctuation approximation. In addition, 
the A-derivative of 5*0 vanishes due to the saddle point condition given in Eq. (4.10). As 
a result, although the saddle point solution (A, ^m) depends on the chemical potential /i, 
one does not need to keep the contribution in Eq. (3.40) in evaluating N — —dQ/djj,, as we 
mentioned in Section III. 

In order to carry out the functional integrals in A in Eq. (4.16), we divide the fluctuation 
flelds into the real part and imaginary part as SA{x) = SAr^x) — iSAj{x) and 6^{x) = 
d^R{x) — iS^i{x). Then we can write SA{x) = A{x)ti + B{x)t2, where A{x) and B{x) are 
real quantities. When we transform SS in Eq. (4.15) from the x{— (r, T))-representation to 
q{= (q, z/„))-representation, we flnd 



- ^E'^J[i + ^n°(g)]x, - lE'Kn^'z)^, + r^n'm. 



9t 
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(4.17) 



Here Xt = {SAR{-q),6Ai{-q)), = (Note that 5A%j{q) = 

5A.R,i{—q) and 5$^j(g) = 5^Rj{—q).) The prime in the ^-summation means that one 
does not take —q when one takes q, in order to avoid double counting. n°(g), W and D^{q) 
are, respectively, defined in Eqs. (3.30), (3.35) and (3.32). Executing the functional integral 
with respect to Xq^ we obtain 

- + - |E'n'n°(g)[i + ^fi°(g)]-^n°(g)f. 



^ETr[l4nO(,)]-|i:'ft 



2/3 



VFZ)°(g)-^i^-^ - ?^fi°(g)[i + ^n°(g) 



-In 



(4.18) 



In obtaining Eq. (4.18), we have used the formula 

^ ' ' det[M] 



(4.19) 



Carrying out the functional integral with respect to Yq using this formula, the second term 
in Eq. (4.16) is given by 



A(A,0^,/.) = ^^Trln[l + |fl°(g) 



+ 4ETrln|-#^°(g)-^l^-^ + |fl°(g) 
-^ETrln[-^^g)-^] + ^ETrln[l + ^[t/-^r'l)°(g)-^]S°(g)]. (4.20) 

The matrix correlation function S°(q, iz/„) is defined in Eq. (3.34). 

Substituting Eqs. (4.13) and (4.20) into Eq. (4.16), the thermodynamic potential Q is 
given by 



1 + |n^,) 



-1 



^^ = ^ + E[ep-^E T>ln[-G'(p, iu;^)-^] 
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(4.21) 



In turn, the total number of particles N — —d^/d/i is given by 

1 ^„ r^. .dG{p)-^ 



-^E|;T.ln[i + >-rf^»M-')=»(, 



-1. 



— tanh 



+ 2l]nB(U) 



- ^ E |;Trln{l + - ^,^D°(q,i^„)]5°(q,z^.„)} 



(4.22) 



This result is identical to that given in Eq. (3.41). The reproduction of the results in Eqs. 
(3.16) and (3.41) [see Eqs. (4.11) and (4.22)] exphcitly proves that the strong-couphng 
theory presented in Section III in terms of summing up diagrams corresponds to keeping 
Gaussian fluctuations around the (mean-field) order parameters A and ^m- 



V. GOLDSTONE MODES IN THE BCS-BEC CROSSOVER REGION 



A. Correlation function and 6-boson Green's function in the HF-RPA 



The Goldstone mode in fermion superfluidity (Anderson-Bogoliubov) is a collective phase 
oscillation (phason) of Cooper-pairs, and thus it appears as a pole in the phase correlation 
function 1122 (q, iJ^n cu + iS). In the present coupled fermion-boson model, we also expect a 
Bogohubov phonon mode associated with the EEC of 6-molecules characterized by the Bose 
order parameter = (fcq=o)- This mode appears in the excitation spectrum of the 6-boson 
Green's function D(q, r) = ~(2^r{^g(''')^g(0)}). However, since the Cooper-pair amplitude A 
and the 6-boson order parameter (p^ are coupled with each other via the Feshbach resonance 
[see Eq. (2.5)], these two Goldstone modes are strongly hybridized. 
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In calculating the Goldstone mode, we have to be careful use a consistent approximation 
for the self-energy and vertex correction. In this regard, apart from the chemical potential, 
we recall that the gap equation in Eq. (3.16) is obtained from the Hartree-Fock Green's 
function in Eq. (3.11). Thus we should employ the HF-RPA formalism for the correlation 
functions ft = {Ilj^} and 6-boson Green's function [24,34,35]. Using the HF-RPA guarantees 
the gapless behavior of the Goldstone mode. 

In calculating the correlation function ft, it is convenient to first consider S = 21^~"^ftW^ 
[see Eq. (3.34)]. Within the HF-RPA, S is described by the sum of the diagrams in Fig. 
5(a). The summation gives 

1 



^ — 



1 

a'- 



1 - -.g'rD'^u 



1 + 



60 



(5.1) 



where S[/ involves the effect of non- resonant interaction —U using the HF-RPA, as shown 
in Fig. 5(c) and given by Eq. (3.36). The correlation function ft is then obtained from the 
inverse relation ft = |l^Sl^~^, namely 



ft = ft° 
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More explicitly, the amplitude and phase correlation functions Hn and II22 are given by 
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Here the lowest-order non-interacting correlation functions 11°^ are defined in Eqs. (3.20)- 
(3.25). Vi describes amplitude- amplitude and phase-phase interactions, while V2 is an 
amplitude-phase coupling mediated by 6-molecules. Since both V2 and vanish at 
q = iun — 0, the denominator of 1122 in Eq. (5.3) reduces to 1 -|- ^1122(0, 0). This vanishes 
when the gap equation in Eq. (3.26) is satisfied below Tc. This proves that the collective 
phase oscillation is always gapless at q = 0, a requirement of the Anderson- Bogoliubov 
(Goldstone) mode. 

The 6-boson Green's function consistent with treating 11 in the HF-RPA is shown in 
terms of diagrams in Fig. 5(b). The result is just the same as the renormalized 6-boson 
Green's function D given in Eq. (3.48). As discussed at the end of Section III, the excitation 
spectrum of D is gapless at q = in the superfiuid phase below T^, which is again consistent 
with the expected gapless Bogoliubov phonon mode. 



B. Goldstone mode 



As discussed in the preceding subsection, the excitation spectrum of the renormalized 
6-boson is determined by 



det[Z)(q, iun^u; + iS)~'^] = 0. 



(5.6) 



On the other hand, the collective phase oscillation and amplitude oscillation are obtained 
from the poles of ft as given by Eq. (5.2) 
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(5.7) 



Comparing Eqs. (5.6) with (5.7), we see that the correlation functions Uij — 1,2) and 
the 6-boson Green's function D have the identical poles (unless det[l + ^n''(q, 0;+)] has ze- 
ros) . This equivalence is due to hybridizing effects from the coupling between the amplitude 
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fluctuations, phase fluctuations and the 6-bosons with the amphtude-phase couphng II12 
and the Feshbach resonance couphng g^.. Thus, in principle, we can consider either one of 
the correlation functions Uij or the 6-boson Green's function D when we want to study the 
Goldstone mode. We work with the phase correlation function 1122 in Eq. (5.3), in which 
case the dispersion relation of the collective mode is given by 

Landau-damping associated with thermally excited fermion quasi-particles leads to the 
imaginary part of collective oscillations and, apart from T = 0, we have to look for a complex 
solution to Eq. (5.8). This requires a complicated analysis. In this paper, for simplicity, we 
only consider the real part of this equation [34,35], 

r Vi~ v;2n?in92 -n°on9, n°2 i , , 

ise[i + f n.. + f "^^^ng - y^T^) = 0. (5.9) 

From this equation, we can obtain real frequencies as approximate solutions. In order to 
check the validity of this prescription, we also solve another approximate equation for the 
mode energy obtained from the renormalized 6-boson Green's function, 

i?e[L>?i(q, iiy 0;+)"^] = 0. (5.10) 

As for the damping of the collective mode, we investigate this effect by examining the 
width of the peak of the collective mode in the spectrum of the correlation functions 

given by /m[nii(q, cl)+)] and /m[n22(q, ^+)], as well as the 6-boson excitation spectrum 
/m[l)ii(q, w„ — > (^+)]. The structure function 

Sjj{ci,uj) = --[nB{uj) + l]Im[fljj{q,iu^uj+)] (j = l,2). (5.11) 

TT 

is more convenient than Im[fljj] in studying collective behavior, because a diffusive mode 
spectrum which arises appears as a central peak at a; = in Sjj, which can be easily 
distinguished from collective modes appearing at a flnite frequency. In contrast, this diffusive 
mode shows up as a peak at a finite frequency ujd in Im[fljj], which exhibits a structure of 
the kind uu / (a;^ + uul) . 
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C. Goldstone mode at zero temperature 



At zero temperature, since Landau-damping is absent for modes below the excitation gap 
2A, we can deal directly with Eq. (5.8) in determining the Goldstone mode for frequencies 
cu below the excitation gap. In the long wavelength or phonon limit, we take cu — v^q and 
expand Eq. (5.8) up to the quadratic order in terms of a; and q. After some calculation, the 
velocity of the Goldstone mode i'^ is given by 
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The factors A and B are obtained from the expansion of 1122 n22(q, Wn) = 1^22(0, 0) + 
Ai/^ + Bq^ /2m, with the explicit expressions 
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Finally, the factor 77 is related to the amplitude-phase coupling 1152, namely 



^ ^ n52(q, i^n) 



T=u„=0 
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(5.14) 



The second term in the denominator in Eq. (5.12) describes the effect of amplitude-phase 
coupling (second order in rj), while the second term in the numerator and the third term in 
the denominator involve the effect of the Feshbach resonance coupling (second order in g,.) . 

1. BCS regime: 2v > 2£f 

In the BCS limit {2v » 2//), the terms involving the factor l/(2i/ — 2//) can be neglected 
in Eq. (5.12). In addition, since the region near the Fermi surface dominates just as in 
ordinary weak-coupling BCS theory, we may take 
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/oo 
dm), (5.15) 
p '°° 

where N{jjL{T — 0)) is the fermion density of states (DOS) at the Fermi surface. In this 

approximation, r] coming from the amphtude-phase correlation function 1152 vanishes, and 

Eq. (5.12) reduces to 

where v-p = y'2m//(0). In evaluating B, we have approximated Ep appearing in Eq. (5.13) 
by the Fermi energy //(O). Eq. (5.16) is the well-known velocity of the Anderson- Bogohubov 
phonon in weak-coupling BCS superfluidity [29,36,45]. 

The same result as in Eq. (5.16) is given by the pole of the phase correlation function 
1122 in the BCS limit. 1122 in this hmit is given by [taking 1152 — 0, Vi — U and = in 
Eq. (5.3)] 

n°,(q,a;+) ^n°2(q,a;+) 

"^^^^'"^^ - TTpiS^ " no,(q,c.^)-no,(o,o)- ^'-''^ 

Using n^2(q)^+) = n^2(0)0) - + Bq^/2m, we easily find that the poles of Eq. (5.17) 
are phonons with the velocity given by Eq. (5.16). This shows that the Goldstone mode in 
the BCS limit is a pure collective phase oscillation of the Cooper-pair order parameter A. 

The amplitude-phase coupling 1152 Eq. (3.8) vanishes in the BCS limit when we use 
Eq. (5.15) and retain up to 0{q/p-p), where p-p is the Fermi momentum. In this case, the 
amphtude mode is decoupled from Eq. (5.8). This mode can then be obtained from the 
amplitude correlation function working within the HF-RPA, 

"^^^^'"-^ - i + fn5,(q,a;,) " n5,(q, c.,) - no,(o, 0) " ^'-^'^ 

In particular, at q = 0, Eq. (5.18) reduces to 

(,^2A). (5.19) 

7rC/Ar(;x(0)) V4A2-u;2 ^ > ^ > 
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Thus, although u — 2A is seen to be a branch cut rather than a pole, a strong peak 
is expected at the excitation gap uj — 2A in the spectrum of the amplitude correlation 
function /m[nii] and hence 5'ii(q, a;) [45,46]. 



2. BEC regime: v < 



In the BEC regime {u < 0), since the chemical potential /i approaches v as ly decreases, 
the effective interaction Ues defined in Eq. (3.2) is dominated by the Feshbach contribution, 
g^/{2iy — 2/x) U. On the other hand, the correlation functions 11°^ become less dominant. 
Since |//| ~ ^ A, the energy gap is also less important in the excitation spectrum of 
fermion quasi-particles, so that we can approximate Ep — [e^ — jiY + A^ ~ £p + In 
this hmiting case, Eq. (5.12) reduces to 
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g'A 



1/2 
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(5.20) 



l32|i/£fP/2- 

where we have rescaled the Feshbach coupling as g^y/N — > g'^ in the last expression. In 
evaluating the p-summation in Eq. (5.20), we have taken into account the correct energy 
dependence of the DOS (oc ^Je^) in contrast to the approximation in Eq. (5.15). Since 
0ni is the dominant contribution to A in the BEC regime, Eq. (5.20) is proportional to 
Y^'iV^ (where — 0^ is the number of condensed 6-molecules) . This dependence on is 
characteristic of the Bogoliubov phonon mode in a Bose-condensed gas. Thus, Eq. (5.20) 
may be regarded as the velocity of the Bogoliubov phonon associated with a condensate of 
6-molecules. We also note that the expression in Eq. (5.20) gives as ~ l/|i/p/^ and this 
approaches zero when v — oo. In this limit, the superfluidity is described by a BEC of a 
free gas with N/2 bosons having the particle-like excitation spectrum uj — g^/2M, with no 
linear (or phonon) component. 

The Bogoliubov phonon associated with a Bose condensate of paired fcrmions has also 
been discussed in strong-coupling superconductivity [44] , as well as for excitons in optically- 
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excited semiconductors [37]. Eq. (5.12) reproduces the Bogoliubov phonon velocity in 
strong-coupling superconductivity [44]. In this case, the Feshbach coupling g,. and the h- 
boson are absent, while the non-resonant interaction —U is taken to be strong. Then Eq. 
(5.12) reduces to 
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(5.21) 



where we have taken the strong-coupling limit {U oo), with // large and negative. The 
last expression can be shown to be equivalent to the Bogoliubov phonon velocity in strong- 
coupling superconductivity, as discussed in Ref . [44] . Using the strong-coupling expressions, 
A = ( 1 6/37r) ^/^£f/ y^ppOs and ji — —l/2ma1 [44], where is an s-wave scattering length of 
the pairing interaction between fermions and p-p is the Fermi momentum, Eq. (5.21) reduces 
to a more familiar form. 



V4> 



M 
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(5.22) 



Here tib = N/2V is the number density of Cooper-pair bosons (where V is the volume of a 
system) and Ub = 47raB/M (M = 2m) is an effective interaction between Cooper-pairs with 
an s-wave scattering length ob = 20^. 



VI. BCS-BEC CROSSOVER IN THE SUPERFLUID PHASE 



In this Section, we present numerical results for the various thermodynamic parameters 
and correlation functions discussed in earlier Sections, as one passes through the BCS-BEC 
crossover region. We give results as a function of the temperature in the superfluid phase 
(T < Tc) for different values of the 6-molecule threshold 2i/. We take the Fermi energy ep, 
Fermi momentum pp and Fermi velocity vp as the units of energy, momentum and velocity, 
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respectively, where £f, Pf, and vp are all in the absence of coupling to 6-bosons. As the 
unit for the number for particles, we take the total number of Fermi atoms A'^. Since the 
interactions U and always appear with N as UN and g^N, we rescale them as UN —>■ U 
and grVN gr- As for the energy cutoff which is necessary in the gap equation in Eq. 
(3.16) and the correlation functions n°^- in Eqs. (3.20)-(3.25), we employ the Gaussian cutoff 
g-Cep/^'^c)^ with ujc = 2£f, as in our previous work for [12-14]. Since our "strong-coupling" 
theory is based on perturbative expansions with respect to —U and g^r, these coupling terms 
are assumed to be weak perturbations. For this reason, we take g^ = O.dep and U — O.Sep 
in all our numerical calculations. In Ref . [14] , we have discussed the BCS-BEC crossover at 
Tc for the case of a very broad [11,15] Feshbach resonance {gr » £f)- 

A. Temperature dependence of the order parameter and chemical potential 

The self-consistent solutions (A,//) of the coupled equations (3.16) and (3.41) are shown 
in Figs. 6 and 7. In Fig. 6(a), we find that the temperature dependence of the order 
parameter A = A — gr4>jn agrees well with the BCS theory ('BCS' in the figure). Since A 
determines the energy gap of the fermion quasi-particle spectrum as Ep — + A^, the 
superfluid character at i'/sf = 2 is found to be the BCS-type. However, Fig. 6(a) also 
shows a sizable difference between A and the Cooper-pair amplitude A — J2p{c_piCp^) ■ This 
is because the BEC order parameter (f)^ — (6g=o) is always induced below Tc due to the 
Feshbach coupling effect in Eq. (2.5), which contributes to A = A — (yfr^m- As illustration, 
when we substitute v/ef — 2, gr/£F — 0.6, U/sp — 0.3 and ji/ep ~ 0.84 into Eq. (2.5), the 
contribution of the condensed 6-boson to A is found to be — 5'r0m = 0.517A at T = 0. This 
means about one third of the order parameter A is associated with (p^ even in the BCS 
regime at v/ef — 2. 

As the threshold energy 2i/ is lowered, the Cooper-pair component A becomes less dom- 
inant while the condensed 6-boson component increases (Fig. 6(a) — > (d)). At the same 
time, the temperature dependence of A deviates from the weak-coupling BCS theory. At 
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v/e-p — —2 in Fig. 6(d), A is dominated by (although the Cooper-pair component also 
exists due to the Feshbach couphng effect). In this case, is well described by the BEC 
order parameter of N/2 atoms in a free Bose gas, given by 

We denote cp^'^ in Fig. 6 as 'BEC 

In the intermediate region of the BCS-BEC crossover, A (and also A and ^m) becomes 
double- valued near T^. For example, in Fig. 6(b), the two self-consistent values for A at 
Tc = 0.213Tf (= TJ", where L stands for lower transition temperature) are zero and 0.14£f- 
In this case, when one decreases the temperature, A jumps abruptly at T^". This is because 
the phase transition is suppressed by superfiuid particle-particle fluctuations, once the phase 
transition occurs, the opening up of the fermion quasi-particle excitation gap (2A) strongly 
suppresses these fluctuation effects, which accelerates the increase of A. On the other hand, 
when we raise the temperature from below, since the superfiuid fluctuations are suppressed 
by the excitation gap, we can exceed T^p staying in the superfiuid phase up to a higher 
temperature (= T^). In Fig. 6(b), we see that = 0.2157^, shghtly higher than T^. At 
T^, A vanishes discontinuously. 

Fig. 7 shows that the chemical potential is decreased as the system approaches the 
BEC regime {v <Q). The temperature dependence of ji is found to be weak below (except 
just below Tc in the crossover regime) compared with the normal phase. As discussed in 
Section III, the chemical potential is temperature-independent below in both the BCS 
limit (// = e-p) and the BEC limit {/i — u). Fig. 7 shows that this feature also holds in the 
intermediate region of the BCS-BEC crossover except just below at least for the model 
parameters we have chosen. 

Now that we have calculated self-consistently the values of the composite order parameter 
A as function of both T and 2z/, we can use the results to discuss the spectrum of the BCS- 
Bogoliubov singlc-particlc Green's function. The excitation spectrum of the BCS-Bogoliubov 
quasi-particles is given by 
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Pf{ui) = -- V/m[Gii(p,ia;^ ^u + iS)], (6.2) 

where Gii{p,u; + iS) is the (ll)-component of Eq. (3.11). When the chemical potential // is 
positive, the excitation spectrum in Eq. (6.2) has a finite gap A as in BCS theory, 



PF(a;,/i>0) = ^ [ ( , , +l)^i2+^u^-A^ 



8£p 



+ ( , - l) V/^ - V'^' - A2e(VA2 + _ ^)le(^ _ A), (6.3) 

where G(a;) is the step function. On the other hand, when jj, is negative in the BEG regime, 
one finds. 



Pf{uj, p<Q)^^ { / ~ + l) V/^ + + ^'©('^ - + /^')- (6-4) 

Va;2 - 



In this case, the Fermi quasiparticle excitation gap is found to be ^A^ + /x^, rather than 
A. This reflects the fact that the threshold energy of free fermion excitations is given by |//| 
when p < 0. In the EEC limit (/^ ~ i/ <^ — ^f), we have 3> A and the excitation gap 
reduces to |//|. 

Fig. 8(a) shows the excitation spectrum of BCS-Bogohubov quasi-particles in the BCS- 
BEC crossover. In the BCS regime > Ef), we find a peak at the excitation edge uj = A 
in the spectrum. This is the well-known coherence peak discussed in the superconductivity 
literature [36], and the quasi-particle spectrum is found to be the BCS- type in this regime. 
This coherence is absent when z/ < 0, where the excitation gap gradually changes from A to 

as the threshold energy 2z/ is lowered. Since ~ i/ in the BEC regime, the energy gap 
+ becomes larger for lower values of v in the BEC regime. 

Fig. 8(b) shows the momentum distribution function of Fermi atoms at T = [36], 
which is given by 

vl^{clc,^) = lil-^). (6.5) 



Since the energy gap A in i?p = \JQ + A"^ is larger for smaller values of the threshold energy 
2i/, the steep decrease of around p{T — 0) &X v/e-p — 2 gradually disappears as 2i/ is 

38 



lowered. In addition, the magnitude of Vp decreases as one approaches the BEC regime 
due to the decrease of the chemical potential /x(T — 0), which reflects the fact that that 
most Fermi atoms form 6-bosons in the BEC regime. The quantity Vp also enters into the 
Bogohubov-transformation [36] to BCS quasiparticles 

ll^ ^ '^pcl^ + VpC-pi, (6.6) 

where is a creation operator of a BCS-Bogoliubov quasi-particle and u^= 1 — v^. 

B. Velocity of the Goldstone phonon mode 

Fig. 9 shows the velocity of the Goldstone mode v ^ at T = as obtained from Eq. 
(5.12). In the BCS regime {u ^ Sp), the mode velocity agrees with the well-known 
Anderson-Bogohubov phonon velocity = vf/V^ in Eq. (5.16). [At u/ef — 2, we obtain 
vf = ^2mfi{0) = 0.92vf for /i(0) ~ 0.84£:f, which gives = 0.53vf.] As the threshold 
energy 21/ is lowered, v,^ decreases sharply and approaches the Bogoliubov phonon mode 
given by Eq. (5.20). Fig. 9 indicates that is strongly dependent on the threshold energy 
2i/ in a uniform Fermi gas. 

Fig. 10 shows the velocity of the Goldstone mode at finite temperatures, obtained from 
Eq. (5.9). In Fig. 10(a), is found to be well described by the Anderson-Bogoliubov mode 
in the weak-couphng BCS theory ('BCS' in the figure) in the whole temperature region, as 
expected for the value u = 2£f- On the other hand, v,^ becomes less than the BCS result 
for the Anderson-Bogohubov mode as the threshold energy 2u is decreased. This is shown 
in Fig. 10(b). Since the order parameter A vanishes discontinuously due to the fluctuation 
effect discussed in the previous subsection, shows a finite jump at Tc in Fig. 10(b). 

C. Dispersion relation of the Goldstone mode 

Fig. 11 shows the dispersion of the Goldstone mode at T = 0.5Tc. In the BCS regime 
(Fig. 11(a)), the gapless dispersion is convex and is confined below the excitation gap at 2A. 
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This convex dispersion relation gradually changes to a concave one as one goes from Figs. 
11(a) to 11(d). In the BEC hmit {v < 0), the dispersion is particle-like (g^), characteristic 
of free 6-bosons. Indeed the dispersion in Fig. 11(d) is well described by a; = g^/2M, except 
for the region of linear dispersion cu — Vtj,q around q = 0. 

In the BCS regime, Eq. (5.10) also has a high-energy solution at to/e-p ~ 2.4, as shown 
in Fig. 12. This solution is also obtained at Tc, as shown in Fig. 13(a). Since the energy 
of this solution is close to the threshold energy of the excitation spectrum of free 6-bosons 
as ^Sg=o — '2v — 2jji — 2.32£f (// ~ 0.84£f), this high energy solution is interpreted as an 
excitation of a 6-boson. 

The high-energy mode solution at shown in Fig. 13(a) may be useful as a guide to 
see how the Goldstone mode changes in the BCS-BEC crossover. Fig. 13(a) shows that the 
energy of this mode gradually decreases as the threshold energy 2i/ is lowered. At the same 
time, other solutions appear around a; = q = 0, as shown for y/e-p — 1.29 in panel (b). At 
v/e-p — 1.14 and below, two finite energy and one gapless excitations combine to produce a 
single phonon. Then the Goldstone mode dispersion below T^. starts to change from convex 
to concave. In Fig. 11(b), although the dispersion below 2 A is still convex, characteristic of 
the Anderson-Bogohubov mode in BCS superfluidity [34] , the dispersion above 2 A obtained 
from Eq. (5.10) approaches the concave dispersion obtained at Tc. In Figs. 11(c) and (d), 
the region where linear dispersion occurs {uj — v^q) is very narrow. The overall behavior is 
the same as the dispersion relation obtained at which is particle-hke, /2M. 

Finally we briefly comment on the approximate Eqs. (5.9) and (5.10). Since the ex- 
citation gap 2 A strongly suppresses the fermion quasi-particle excitations far below Tc, 
the intraband terms of 11°^ and also the Landau damping below 2A are less important at 
T = 0.5Tc. As a result, Eqs. (5.9) and (5.10) both give good approximations to the solution 
of Eq. (5.7) at T = 0.5Tc, and give almost the same results below 2A, as shown in Figs. 
11(a) and (b). On the other hand, the interband terms given by the second hues in Eqs. 
(3.20)-(3.25) give rise to Landau damping when uj > 2A. Thus, above 2A, Eqs. (5.9) and 
(5.10) may have different solutions, because they do not include the imaginary part in the 
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same way. Indeed, the dashed hne above 2 A shown in Fig. 11(b), which is a solution of Eq. 
(5.10), is not obtained from Eq. (5.9). In such a case, a more careful analysis is necessary to 
obtain the correct dispersion relation. However, since the fermion quasi-particles are absent 
in the BEC hmit {u <0), the Landau damping becomes weak in the BEC regime even near 
Tc. Thus Eqs. (5.9) and (5.10) again give almost identical results in the whole energy region 
shown in Figs. 11(c) and (d). 

D. Spectral weight and damping of the Goldstone mode 

In Sections VLB and VI.C, we considered the Goldstone mode neglecting the Landau 
damping. In this section, we evaluate the damping from the width of the collective mode in 
the structure function Sjj{q,u!^) in Eq. (5.11) as well as in the 6-boson spectral density 

Pb{(1,uj) = --Im[bn{c[,uj+)]. (6.7) 

TT 

In Figs. 14-17, we show the 6-boson excitation spectrum, as well as the phase {S22) 
and amplitude (Su) structure functions for T, below Tc. In Fig. 14, we find that the 
Anderson-Bogoliubov Goldstone mode does not appear as a visible peak in the spectrum 
at T/Tc = 0.99. (The inset in Fig. 14 shows that the mode energy is uj/sf = 0.01 for 
q = 0.02pF-) This means that the Anderson-Bogoliubov mode is over-damped near Tc 
because of strong Landau damping from the thermally-excited fermion quasi-particles. We 
note that /S'22(q, uj) and 5*11 (q, uj) both exhibit strong central peaks at a; = for T/Tc = 0.99, 
which indicates the presence of a large number of thermally-excited fermion quasi-particles 
expected at this temperature. 

At low temperatures, when the Landau damping becomes weaker, the Anderson- 
Bogoliubov mode appears as a visible peak in ps{q,u!) and S22{(l,uj) as shown in Fig. 
14. The peak width becomes narrower at lower temperatures, reflecting the weaker Landau 
damping by fermions. (At T = 0, it becomes a sharp 5-function peak.) The peak position 
{uj/e-p ~ 0.021) at T/Tc = 0.5 agrees well with the dispersion relation shown in Fig. 11(a). 
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Since the Anderson-Bogoliubov mode is a collective phase oscillation of the Cooper-pair 
order parameter A, the appearance of this mode in pB{q,uj) indicates the presence of the 
coupling between the 6-bosons and the phase fluctuations of A. On the other hand, no 
peak structure is observed except for the central peak at a; = in the amplitude structure 
function Su shown in Fig. 14. Only a shght structure appears at oj/ef ~ 0.02 at T/T^ — 0.7, 
as shown in part (a) of Fig. 18. (The peak at lu/ef — 0.07 in Fig. 18 is the amplitude 
mode at the edge of interband excitations u; — 2 A.) This is because the amphtude-phase 
couphng is very weak in the BCS regime, so that the collective phase oscillation does 
not strongly couple into the amplitude fluctuations. 

As shown in Fig. 12, Eq. (5.10) also has a high energy solution in the BCS regime at 
u/ef — 2. This solution is the strong resonance in pB(q, '^), as shown in part (b) of Fig. 18. 

In Fig. 15, we find a broad peak in 5*22 at T/T^. = 0.99. As shown in the inset in Fig. 15, 
the peak position is different from the one expected from Eqs. (5.9) and (5.10). However, 
since this resonance shows gapless behavior, it clearly must be the Goldstone mode. (The 
difference shown in the inset is due to the Landau damping effect.) Indeed, the peak energy 
at T/Tc — 0.5 agrees well with the dispersion in Fig. 11(b). As expected. Landau damping 
from fermions becomes weak as one approaches the BEC regime. When v < 0, we can 
observe the sharp peak structure even near Tc in Figs. 16 and 17, and the peak position 
always agrees well with the dispersion relation given in Fig. 1. The dispersion relation of 
the Goldstone mode approaches the temperature-independent particle-like one in the BEC 
regime, and hence the temperature dependence of the peak energy in Fig. 17 is weak. 

We also find in Figs. 15-17 that the Goldstone mode appears in the amplitude structure 
factor 5'ii(q, a;), which indicates that the amplitude-phase coupling becomes stronger as 2i/ 
is decreased. Although the amplitude fluctuations described by IIn are not important in 
the BCS regime, we cannot neglect these fluctuations in the BCS-BEC crossover regime. 



42 



VII. DENSITY-DENSITY CORRELATION FUNCTION 



An important problem is how to experimentally observe the Goldstone mode discussed 
in Section V and VI. In this section, we show that the density-density correlation function 
1133 defined in Eq. (3.19) exhibits this mode as a pole. This correlation function can be 
experimentally probed by many techniques, including two-photon Bragg scattering [50]. 

Since the Goldstone mode is associated with the collective phase fluctuations of the 
Cooper-pairs, a coupling between density fluctuations and phase fluctuations is necessary in 
order for the Goldstone mode to appear in the spectrum exhibited by Il3^{c[,uj+). In this 
regard, we note that a phase-density coupling exists in fermion superfluidity because of the 
presence of the Josephson effect. This couphng is described by 1123 (and 1132) defined in Eq. 
(3.24). In contrast to the amplitude-phase correlation function n5'2, which is very weak in 
the BCS regime, 1123 is finite even if one works with the approximation in Eq. (5.15). Besides 
1123, density fluctuations also couple with superfluid fluctuations through amplitude-density 
coupling n'j'3 (and H^i) defined in Eq. (3.25). This coupling is not important in the BCS 
regime, because 11^3 vanishes when one uses the approximation in Eq. (5.15) and neglects 
term of order 0{{q/pFY)). However, as in the case of the amplitude-phase correlation 
function 1152, ^^is coupling effect becomes stronger as one goes into the BCS-BEC crossover 
regime. Then the density-density correlation function 1133 is found to couple with superfluid 
fluctuations (amplitude and phase) via 1123 and 1153, as shown diagrammatically in Fig. 19. 

The density-density correlation function 1133 given by the HF-RPA is obtained by extend- 
ing the method discussed in Section V [40]. When we introduce a 3 x 3- matrix correlation 
function II = {11^} — 1, 2, 3), which involves the density-density correlation function as 
the (33)-component, we obtain an equation similar to Eq. (5.2), namely 



n = flHl 








,1-1 



(7.1) 



Here the interaction 3 x 3-matrix U is defined by 
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^ f/ll Ul2 ^ 



u 



U21 U22 




(7.2) 



where Uij = {U — g^WD^W~^}ij. If we included a density-density interaction term to our 
coupled fermion-boson model in Eq. (2.1), it would be included as the (33) -component in 
Eq. (7.2). The density-density correlation function is obtained from the (33) -component of 
Eq. (7.1), 



n 



33 



n^3 + 



i + f[n?, + no 



22J 
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n^i] + fr + ^][n?ino 



22 



n?2n^i] 



(7.3) 



where the numerator is 



n33 = + - n°,n°,n?3 + ul,ul,u% - n°,n2,n: 



23] 



yin^inS3 



■'-'■32-'-'-23J 



y[n^in^3 



(7.4) 



Each term in Eq. (7.4) involves correlation functions between density fluctuations and 
superfluid fluctuations, such as II^s and Ul^. 

Fig. 20 shows the dynamic structure function related to the spectral density of density- 
density correlation function. 



S'33(q,a;) = ~-[nB{ui) + l]/m[n33(q, ^ uj+)]. 

TT 



(7.5) 



At low temperatures in Figs. 20(a) and (b) (BCS regime), and at all the temperatures in 
Figs. 20(c) and (d) (BEC regime), we can clearly see the Goldstone mode as a resonant peak 
in 5'33(q, a;). The peak position is the same as that in the spectra shown in Figs. 14-17. 
Since the density-density correlation function 1133 can be measured more easily than the 
phase correlation function 1122, and may be the most useful way of observing the Goldstone 
mode. 
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VIII. SUMMARY AND DISCUSSION 



In this paper, we have investigated the BCS-BEC crossover in the superfluid phase of a 
uniform gas of Fermi atoms with a Feshbach resonance. Going past the simple weak-couphng 
mean-field theory, we included the strong-coupling effect originating from the pairing interac- 
tion associated with the Feshbach resonance. We have extended our previous work [12-14] at 
and above to the superfluid region below Tc. We showed that the superfluid order param- 
eter continuously changes from the Cooper-pair amplitude A = U J2p{c-piCpi) in the BCS 
regime to the square root of the number of condensed 6- molecules = (6,=o) associated 
with the Feshbach resonance in the EEC regime, as one lowers the threshold energy 2i/ of the 
Feshbach resonance. In the intermediate regime in the BCS-BEC crossover, superfluidity is 
described by the composite order parameter A = A — gr(l)m [8,9,12]. 

The Goldstone mode is one of the most fundamental phenomena in an ordered system 
with spontaneous breakdown of a continuous symmetry. In this paper, we investigated how 
the Anderson- Bogoliubov mode, which is the Goldstone mode in BCS superfluidity, changes 
to the Bogoliubov phonon mode in the BCS-BEC crossover. The velocity of the Goldstone 
mode Vfj) strongly depends on the threshold energy 2i/, and decreases as one approaches the 
BEG regime. The Anderson-Bogoliubov phonon mode may be a useful way of monitoring 
the BCS-BEC crossover. Since it is difficult to strongly modify the strength of the pairing 
interaction in metallic superconductors, the tunable pairing interaction associated with the 
Feshbach resonance in Fermi atomic gases gives one a unique tool to clarify the physics in 
the BCS-BEC crossover region. 

We also investigated the damping of the Goldstone mode. The fermion Landau damping 
of the Goldstone mode due to couphng to fermions becomes weak far below Tc, reflecting the 
fact that the thermal excitation of fermion quasi-particles is negligible. This damping effect 
is always weak in the BEG regime because the system is then dominated by 6-molecules, 
with the suppression of the Fermi quasiparticle spectrum. Thus, except in a small region 
near Tc in the BCS regime, the Goldstone mode appears as a strong resonance both in the 
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spectrum of the phase correlation function and in the excitation spectrum of 6-bosons. As 
a way to experimentally observe this Goldstone mode, we noted that since the amplitude 
and phase fluctuations couple with the density fluctuations, one can observe this collective 
mode in the spectrum of the density-density correlation function 1133 (q, In cigar-shaped 
trapped Fermi gases, the direct observation of a density fluctuation pulse might be possible, 
in analogy with the observation of a Bogoliubov phonon in a superfluid Bose gas with a 
very weak axial trapping potential [47]. We also note that two photon Bragg scattering 
experiments provide a convenient way of studying density fluctuations in a trapped atomic 
gas [50,51]. In such experiments, /mn33(q, a;) is measured directly, rather than 5'33(q, a;) 
[52]. 

It is important to remember that our treatment of the BCS-BEC crossover leaves out 
several important contributions. In our calculation of various response functions, we always 
ignored the fermion self-energy arising from the coupling to the non-Bose-condensed bosons 
[38] . In the crossover region, these two-particle states are strongly damped. The numerical 
results given in Section VI and VII show how Landau damping due to fermions decreases 
as we approach the BEG regime (small or negative values of vje-p). This is simply because 
more and more the fermions are forming bound states. However, we expect new forms of 
damping to arise in the BEG region, namely the Landau and Baliaev damping associated 
with the interaction between Bogoliubov excitations [48]. 

In this paper, we have not considered the effect of a trap. The atomic density profile in 
the BCS-BEG crossover in a trap was recently investigated within the LDA by the authors 
at and above [13,14], using the coupled fermion-boson model, and at T = by Perah et 
al. [28], using the strong-coupling BGS model [26]. A trap potential also leads to various 
discrete low energy collective modes associated with the confined geometry, and it will be 
interesting to study how these collective oscillations behave in the BGS-BEG crossover. This 
will be the subject of a future paper, but we conclude with a few brief remarks on this. In the 
BGS limit of our model, one has the extensive theoretical work [31,32,48] on the Gooper-pair 
condensate modes in trapped Fermi gases. To be specific, we consider the quadrupole mode 
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cuq — y/2u!o (where cuq is the trap frequency of a spherical trap) . This is an analogue of the 
Goldstone phonon mode in a uniform system. If ujq is greater than the effective threshold 
2Aefr for breaking up Cooper-pairs, this mode is damped and has small spectral weight [32]. 
However, as one approaches the BCS-BEC crossover, the spectral weight of the fermion 
quasi-particles decreases, shifting to the quadrupole mode cuq. In the BEC limit, the mode 
cuq = \/2u!o is the well known quadrupole oscillation of a trapped Bose gas [48,49]. In this 
paper, we have seen that the Goldstone phonon mode in a uniform gas also persists in the 
BCS-BEC crossover, but the phonon velocity changes. In contrast, the quadrupole mode 
frequency does not depend explicitly on the interaction strength, having the same frequency 
loq — \/2ujq in both the BCS and BEC hmits. What does change is its spectral weight and 
damping, which could be used as an experimental signature. 
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FIGURES 

FIG. 1. (a) t-matrix approximation for the particle-particle scattering vertex F at Tc. The solid 
line represents the fermion Green's function. The first line in the figure involves the ladder processes 
by non-resonant interaction —U, while the second line involves the Feshbach resonance described 
by the 6-boson Green's function Dq. (b) Fluctuation contribution to the thermodynamic potential 
O at Tc. The first line represents the Cooper-channel particle-particle fluctuations associated with 
the non-resonant interaction — f/, and the second line describes the effect of the Feshbach resonance 
coupling Qt-. (c) The shaded bubble includes multi-ladder scattering processes by —U. 

FIG. 2. (a) BCS-BEC crossover at the superfiuid phase transition temperature in a gas of 

Fermi atoms with a Feshbach resonance, as a function of the threshold parameter v. We take 
gr/^F = 0.6, U/ep = 0.6 and Uc/e-p = 2. 'BCS' labels the result in the absence of fluctuation 
effects, and 'EEC' gives of a Bose-condensed gas of iV/2 molecules of mass M = 2m. (b) The 
chemical potential ji at Tc, shown as a function of v. These results are from Ref.l2. 

FIG. 3. The off-diagonal static mean-field self-energies in the 2 x 2-matrix fermion Green's 
function G(p, ia;^) given in Eq. (3.11). Diagram (a) gives the contribution from the non-resonant 
interaction —U. Diagrams (b) and (c) include the pairing interaction mediated by a Feshbach 
molecular boson described by the 6-boson Green's function Do(q, ^t'n)- In the diagrams (b) and 
(c), r± = Ti ± iT2- Since G(p,iu;) in Eq. (3.11) does not have a r2-component, a diagram similar 
to (a), where ri is replaced by r2, is absent. 

FIG. 4. (a) Fluctuation contribution ((50) to the thermodynamic potential below T^. The 
bubble shows the correlation function 11?^. (b) Correlation function of amplitude fluctuations 
IIii involving a coupling with phase fluctuations described by the response function 1122- The 
shaded bubble includes multi-scattering processes by —U. Since the 2 x 2-matrix fermion Green's 
functions G are given in the Nambu representation, 11°- is formally described as a particle-hole 
bubble diagram. (In contrast, 11 is described by a particle-particle bubble diagram in Fig. 1 
because the Nambu two-component representation is not used there.) 
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FIG. 5. The correlation function E!(q, zz/^) is shown in (a) and 6-boson Green's function 
i)(q, is shown in (b), both within the HF-RPA. The shaded bubble S(7(q, shown in 
(c) includes RPA-type diagrams involving the non-resonant attractive interaction —U. 

FIG. 6. Temperature dependence of the superfluid order parameter A = A — 5r<^m in the 
BCS-BEC crossover region. We also show the separate Cooper-pair component A and the 
Bose-condensed 6-molecule component (^m- We take gr/sp = 0.6 and U/ef = 0.3, while A and A 
are normalized with respect to ep. The character of the superfluidity changes from the BCS-type 
to the BEC-type as one goes from (a) to (d). In this figure, 'BCS' labels the order parameter given 
by the weak-coupling BCS theory which omits the particle-particle fluctuations. 'BEC labels the 
order parameter (j)^ of a free Bose gas of N/2 atoms, as given by Eq. (6.1). 

FIG. 7. The chemical potential /Lt as a function of temperature, for gr/sp = 0.6 and U/ep = 0.3. 
The solid circles show the superfluid phase transition temperature. 

FIG. 8. (a) Density of states of BCS-Bogoliubov quasi-particles at T = 0. For u/sp > 1, 
the threshold energy of the quasiparticle excitations is given by w = A, at which a coher- 
ence peak appears. For v/ep < 0, the coherence peak is absent and the spectrum starts from 
cj = (/x^ -I- A^)-*^/^. The excitation density of states approaches that for a free Fermi gas of atoms 
pp{uj) = (3 A^/4ep''^ ) + /i)^/^ in the high-energy region, (b) Momentum distribution function 
Vp = {c^pcrCpa) of Fermi atoms with a-spin component at T = 0. The solid circle indicates //(T = 0). 

FIG. 9. Velocity of the Goldstone phonon (normalized to the Fermi velocity) at T = in 
the BCS-BEC crossover region. 

FIG. 10. Temperature dependence of the Goldstone phonon velocity v^, (a) v/ep = 2 (b) 
u/ep = 1. 'BCS' labels the Anderson-Bogoliubov phonon velocity given by the weak-coupling BCS 
theory, where we use Tc from Fig. 2 and fj, from Fig. 7. 
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FIG. 11. Dispersion relation of the Goldstone mode at T = O.STc, for different values of v. 
In this figure, the curves labeled as [Tc] shows the dispersion relation at obtained from Eq. 
(5.10). In panels (a) and (b), Eqs.(5.9) and (5.10) give almost the same results below u} = 2A. At 
v/e-p = 2, Eq. (5.10) has another high energy solution around u/e-p = 2.4, as shown in Fig. 12. 

FIG. 12. The high energy solution of Eq. (5.10) in the case of v/sp = 2, for T = Q.bT^. The 
lower energy solution is the Goldstone mode shown in Fig. 11(a). 

FIG. 13. (a) Excitation spectrum of the renormalized 6-boson Green's function at Tc obtained 
from Eq. (5.10). At u/e-p = 2, there is another solution at w = q = corresponding to the 
Goldstone mode. At v/ep = 1.5 and 1.29, two finite energy modes and one gapless mode are 
present, two of them near q = a; = 0. For illustration, we show these low energy modes for 
u/ep = 1.29 in panel (b), where the two low energy solutions come together and disappear at 
q/pp ~ 0.074. At v/ep = and 1.14 in panel (a), only one gapless mode is present. 

FIG. 14. Spectral weight of the renormalized 6-boson, pB(q, ^^) = — ^-f'Ti[£)ii(q, w)], and the 
structure functions S'ii(q,a;) and (q,w). These results are for v = 2£p (in the BCS region) and 
q = 0.02;jF- The inset shows the dispersion relation of the Goldstone mode at T/Tq = 0.99: (1) 
and (2) are obtained from Eqs. (5.9) and (5.10), respectively. 

FIG. 15. Same plots as in Fig. 14, for u = £p and q = 0.15^?f- In the inset, the 'peak energy' 
gives the peak position in S'22(q, w) at T/Tc = 0.99. 

FIG. 16. Same plots as in Fig. 14, for z/ = (crossover region) and q = 0.3pp. 

FIG. 17. Same plots as in Fig. 14, for v = —ep (BEC region) and q = 0.3pp. 
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FIG. 18. The amplitude structure function S'ii(q, w) is shown in (a) and the 6-boson spectrum 
PB(q, w) is shown in (b), for u/ey = 2 (BCS region) and T/Tc = 0.7. In panel (a), the extremely 
small structure visible at w/ep — 0.02 is due to the Anderson-Bogoliubov mode. The larger peak at 
u/e-p — 0.07 is the Amplitude mode. However, since the excitation gap is also at 2 A = 0.068, this 
mode coincides with the edge of the interband excitations. In panel (b), the peak on the left is the 
Anderson-Bogoliubov mode, while the broad peak at w/ep corresponds to the high energy solution 
shown in Fig. 12. A very small peak at uj/e-p ~ 0.068 in panel (b) is located at the excitation gap 
2A. 

FIG. 19. The density-density correlation function 1133 coupled with superfluid amplitude 
and phase fluctuations described by 11 jj = 1,2) (shaded bubbles). 1123 and Ilgg describe 
phase-density correlations, while 1153 and li^^^ give the amplitude-density correlations. 

FIG. 20. Spectrum of the dynamic structure function >S'33(q, w) in the BCS-BEC crossover. 
The momentum values are the same as in Figs. 14-17. The sharp peak structure is the Goldstone 
phonon mode. 
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